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ON THE FLUTTER OF CYLINDRICAL SHELLS AND 
PANELS MOVING IN A FLOW OF GAS* 

By R . D. Stepanov 


The equations of shells are taken in the form, of the general techni- 
cal theory of shallow shells (ref. l) and shells of medium length of 
V. Z. Vlasov (ref. 2). The aerodynamic forces acting on a shell are 
taken into account only as forces of excess pressure according to the 
formula proposed hy A. A. Iliushin (ref. 3 ) • In this work the following ' 
notation is used: a and p are dimensionless coordinates of a point 

on the cylindrical surface of the shell. The coordinate a,, expressed 
in terms of radius R, represents the distance along a generator, and 
p represents the central angle. The dimensions R, h, and Z are 
the radius, thickness, and length of the cylindrical shell; E, a, and 
p are Young's modulus, Poisson's coefficient, and the density of the 
material of the shell, and D is the cylindrical rigidity. The quanti- 
ties u, v, and w are components of the vector of displacement of the 
shell, V is the velocity of flow, v Q is the velocity of sound at 

infinity, p Q is the pressure of the gaseous medium at infinity, and k 

is the exponent of polytropy . The symbol co = p + iq is the complex 
frequency, and c and B are constants . Also, 



D = 


Eh^ 

12(1 - a 2 ) 



da 2 dp 2 


The quantity 


is the coefficient of damping, and Z is the trans 


verse component of the load. 


1. INITIAL RELATIONS OF THE THEORY OF CYLINDRICAL SHELLS 


In the system of dimensionless coordinates a, P, in the case when 
at each point a load, directed along the normal to the surface acts on 

*"0 Flattere Tsilindricheskikh Obolocliek i PaneleT, Dviziiushchikhsia 
V Potoke Gaza." Prikladnaia Matematika i Mekhanika, vol. 21, no. 5, 1957, 
PP. 6M1--657. 

NACA Reviewer 1 s note: The original Russian publication contains cer- 

tain typographical errors and obvious omissions in equations that have been 
corrected without comment. 



2 


NACA 5PM 1438 


the shell (X = Y = 0, Z ^ 0.), the fundamental equation of shallow cylin- 
drical shells in a form convenient for solution is (ref. l) : 


V 2 V 2 V 2 V 2 4 + 1 - = *£ Z 

c 2 jj ' 


( 1 . 1 ) 


where <&(a,p) is a scalar function, defined according to the formulas 


u = 


„ d?<b 


da dp 


- a 

2 bet? 


v = - 


£* + ( 2 + „) Jb- 

dp 5 da 2 dp 


w = V 2 V 2 4 


( 1 . 2 ) 


The internal forces of shallow cylindrical shells are defined through* 
the function <t(a,p) by the following group of formulas: 


H _ Eh d h <b 
R da 2 dp 2 

Ml “B^(& + "S) V2v2t 

11 


s __ Eh 

da^dp 

► 

H - - ~ g ) v 2 V 2 $ 

r 2 5a dp 

. "i v 2 V 2 7 2 « 

r 5 da 

Q 2 . . D a v 2 2 2 
r 3 dp 




Generalized transverse forces defined in the sense of KLrchoff being 
necessary for the formulation of the boundary conditions, are computed 
according to the formulas 
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The internal forces of cylindrical shells of medium, length are 
expressed through the function hy the following group of 
formulas : 



h 
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*L = 


~ ll 

Eh H 

TT _ Eh 

. c 2 

R aa 2 8p 2 

iH 2 - R 

1 - a 

D — 




R 2 


Mo = — 

* R 2 I 





Eh ^®1 

W __ JD (■> _ rr\ 

a 6 ®- 

R 3a 5 ap 

0 ' 

R 2 

, W / 

da 8f 

D 


0 __ D_ 

+ 

d\ 

~"r5 

_cfcc 8p^ 8a 8p^_ 

% ■ r5 

_a P T 

ap5_ 


> (1-T) 


In each particular case, it is necessary to adjoin given boundary con- 
ditions on the edges of the shell to the differential equations (l.l) 
or (1.3). 


2 . STATEMENT OF THE PROBLEM 


The expression for the transverse load Z acting on an element of 
the surface of the shell is composed of two parts: the force of inertia 

Z L = -ph & ' (2.1) 


and the force of aerodynamic action of the supersonic flow directed 
along a generator and flowing around the shell on the outside, which 
is taken into account according to the formula proposed by A. A. Iliushin 
(ref. 3) : 


Z 


_ BV 8w -d cto 

2 " T " Bl St 


(2.2) 


Substituting expression (2.2) into equation (l.l) and taking into 
account the third of the relations in the group (l.2) yields the differ- 
ential equation of small vibrations of shallow cylindrical shells: 


c* 2 V 8 G 




td2 

+ T P 


8t 2 


4 - 


BVR 

Eh 


|L ® 


r 2 b x 

Eh 




(2.3) 
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Correspondingly, substituting the expression for the external load Z 
in equation ( 1 . 5 ) and taking into account the third relation of group (l.6) 
gives the equation of small vibrations of cylindrical shells of medium 
length: 


, 

8 a 4 


2L+ II 

8 p 2 


2 


8 6 


n 2 

+ p 

bp 14 - E St 2 ^ 


*1 _ S3E. 


B X R 2 


Eh 


8 a 8^ 


Eh 


8 t 8p 


= 0 


( 2 » 


The nev dimensionless quantity 


c * 2 -— ! £ (2.5) 

1 - cr 2 12R 2 (l - a 2 ) 

is introduced in equations (2.3) and (2.1+) . 

In all of the following calculations B-]_ will be taken equal to B. 

In 1951+, under the guidance of A. A. Iliushin, an investigation was 
made (refs, k and 5 ) of the self-induced vibrations of a plate moving in 
a gas, which defined in many respects an approach and methods of solution 
of the problem set down; certain results borrowed from the indicated 
works will be introduced below without derivation. 

Examined in this paper is a class of solutions of the form 

®(a,p,t) = ^(ajpje^ (2.6) 

where cu = p + iq is a constant complex frequency. 

In the class of solutions ( 2 . 6 ), the problem of flutter consists 
in determining the least velocity of flow (this velocity will be called 
the critical velocity), which on being exceeded would result in a posi- 
tive real part of the complex frequency. 

Substituting equation ( 2 . 6 ) into equation (2.3) gives, after can- 
celling e^*, 


Here 


2^8^ + akfc _ >vvH - — vH = 0 
8a 4 Eb ^ 




( 2 . 7 ) 


( 2 . 8 ) 
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Equation (2.8) permits us, for each particular eigenvalue A, two 
values of frequency to he defined: 



(2.9) 


For certain eigenvalues, let one of the roots of equation (2.9) he 
a pure imaginary number . Then from equation (2.8) it is easy to obtain 


Re A = \ = p j q 2 


1,0 * = ’v --fr « 


( 2 . 10 ) 


On the complex plane equation (2.10) represents the points 

of the square parabola (fig. 2)* 


A 1 = p^A2 2 (2.11) 

B 2 R 2 

which, following the example of references 4 and I s called the para- 
bola of stability. The region lying inside the parabola of stability 
corresponds to the eigenvalues for which the rootB (2.9) have a negative 
real part, "while the region lying outside the parabola corresponds to 
eigenvalues for which the real parts of one of the roots (2.9) has a 
positive real part. 


3. UNBOUNDED CLOSED CYLINDRICAL SHELL 


For the case of an unbounded closed cylindrical shell, the solution 
of the fundamental differential equation of small vibrations (2.7) will 
be sought in the form 


*(»,P) 


00 00 


XI 

n=l k=l 




(5.1) 


-ft. 

Translator’s note: This apparently refers to figure 2 of refer - 

erences 4 and 5 which is included with the figures of this paper for 
the convenience of the reader. 
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where n is a certain constant number and n and k are constant 

numbers denoting the number of half waves in the meridional direction 
and in the direction of the generators of the shell. 


Substituting equation (3*1) into equation ( 2 . 7 ), yields a charac- 
teristic equation from which the following expression for A results: 


A-j_ = c* 2 (k 2 + n 2 ) 2 + 




(k 2 + n 2 ) 


Ag = -g*k 
2 Eh 


(3-2) 


On the complex plane A;j_,Ag equation (3 .2) represents points of a 
parabola of the eighth power: 


Ai = c*2 


■ *2 + n 2 
b 2 r 2 v 2 


E^h^Aa^ 

E2h2 


B^V 2 


Ag 2 + n c 


-2 


(3.3) 


For the determination of the critical velocity of flow, an analysis 
is made of the problem of the mutual arrangement of the parabola (3-3) 
with the parabola of stability ( 2 . 11) in the cases n = 0 and n ^ 0. 
For n = 0 (that is, for the case when the contour of a transverse sec- 
tion of the shell remains a circle in the process of deformation) equa- 
tions (3.2) take the form: 


A]_ = c* 2 ^ + 1 




*2 

^ Eh 


(3 A) 


For the points of mutual intersection of the parabola (3 A) with 
the parabola of stability, the equalities 


p = c * 2 k^ + 1 

E * 

BR 2 , .. BVR ,, 
Eh 4 ” Eh 


(3.5) 


are valid. 
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Eliminating from the first of equations (3*5) the parameter q 
yields one equation for the determination of the points of mutual inter- 
section of the two parabolas being investigated: 


_ 0 


Ec/ 


= 0 


(3.6) 


the solution of which will he 

f 

pV 2 

k l,2,3,4 = ± \ ~' ~ i 1 
’ 2Ec* 2 


( pV 2 > 

l2Ec* 2 ; 


i 1/2 




From equation (3-7), it follows that, for 

/ 2Ec*\l/2 


> 


* \ P 


(3.7) 


(3.8) 


the parabola (3.4), intersecting with the parabola of stability in four 
'points, extends outside the domain of stability. Hence, it follows ‘that 

for a velocity of flow larger than (2Ec Jt /p) 1 / 2 7 the motion of the shell 
must be unstable. . .. — 


For the investigation of the problem of the mutual intersection of 
the parabola of stability with the parabola (3 • 3 ) ^ in the general case 
for n f 0 the following equation is obtained: 


k 8 + k 6 (kn 2 - p -iSL.'l + k 1 * (fin 4 + -i_ - 2 p jM'l + 


\ 


EC. 


V 


Ec. 


- pn^ - -r- \ + n® = 0 


Ec. 


(3.9) 


Solution of equation (3.9) 


r 


kl = ±2n 2 


[-a ± ( 



gives the eight roots: • 


a 2 - 4b + 8n^) 1//2 ] 
8ng l/2 ) 2 - l6n^j * 


- 1/2 



( 3 . 10 ) 
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where 



b 


2n^a + 




(3.11) 


Similarly to that which was done above for the case n = 0, it is 
possible here also to show that the necessary and sufficient conditions 
for which the parabola (3*3) intersects with the parabola of stability, 
and hence falls outside the domain of stability, reduce to the deter- 
mination of the condition of the appearance of complex roots of equa- 
tion (3.10). 


\ 


Analyzing the expression (3-10), it is possible to set down the two 
following essentially distinct necessary and sufficient conditions that 
the parabola (3-3) intersect with the parabola of stability, and hence, 
extend beyond the domain of stability: 

a^ - 4b + 8n^ = - - -it— = 0, -a t (a^ - 4b 

<=. 2 

In expressions (3-12) it is necessary to satisfy 


+ 8nM l/2 < 

( 3 - 12 ) 

the inequality 


a = 


P 


V s 

Ec* 2 


< 0 


(3.13) 


The inequalities (3-12) and (3-13) make it possible to determine 
the critical velocities: 



n < \ (c*)" 1 / 2 = n* 


V 2 ifl(l6n^c* 2 + 
* 2n[p v * 



(3.1*0 

(3.15) 


The formula for the critical velocity (3*1*1-) identically coincides 
with the critical velocity of flow found for a closed cylindrical shell 
with n = 0, and, as is seen from the inequality (3-13), it may be used 
for all values n < n* which, for thin shells, corresponds to the number 
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of half waves n « 30 to 50, that is, to such a large number of half- 
waves that the form of the transverse section differs little from a _ 
circle . 

The minimum velocity (3 .15)* according to n, takes place for 

n = i c*- 1 / 2 and exactly coincides with the critical velocity found 
above for n = 0. 


From this analysis it follows that flutter of a closed cylindrical 
shell of infinite length in a supersonic flow can possibly take place 

for velocities of flow V > (2Ec # /p)^ 2 when the form of the transverse 
section remains circular. 


Using the formulas (2.9) and (3.2), it is possible 
values of frequency whose essential form depends on the 
flow: 


r 





to obtain two 
velocity of 


+ 


BVk 

ph 



(3.16) 


The solution of the differential equation of small vibrations of 
shallow shells (2.3), taken in the form 


®(a,0/t) = e 1 ( n P +ka )e(P +1< l) t (3.17) 

signifies that along the generators of the shell are propagated waves 
traveling with the velocity 


v b 



(3.18) 


Separating the real part of the complex frequency (eq.. (3 .16) ) from 
the imaginary part yields 


v b 



J 


+ 


B^k 2 

h 2 p 2 R 2 [ 



J 


(3.19) 
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Using formula (3 .19) > the velocity of propagation of a traveling 
wave is determined for V = 0: 


4 * 

E 

c 2 + n 2 ) 2 + 

_k 2 _ + B 2 ~L 

v b ~~ 

pR 2 

_* k 2 

(k 2 + n 2 ) 2 4p 2 h 2 k 2 J_ 


1/2 


The minimum velocity of propagation of a traveling wave occurs for 


n 2 = k 


1 \Va 




■ ' - k 


( 3 . 20 ) 


Here, it has the value 


b min 


SBc * 2 


B 2 


_ pR 2 4p 2 h 2 k 2 


1/2 


( 3 - 21 ) 


h. UNBOUNDED CYLINDRICAL PANEL, SIMPLY SUPPORTED ALONG ITS 

GENERATORS 


In the case of an unbounded cylindrical panel, simply supported along 
it6 generators, the differential equations of small vibrations (2.3) must 
be accompanied by the boundary conditions; 

u = w = 0 N 2 = M 2 = 0 at p = 0 and p = = s/R 


Defining, according to formulas (l.2) and 1.3), the displacements 
and internal forces of the shell through a potential function ®, it is 
possible to write the boundary conditions on the edges (3 = 0, p = p^ 

in the form 


a 2 ® = sS _ B 6 ® _ 
2 Bp 4 ~ 5j? ~ 


Representing the solution of the indicated boundary-value problem 
in the form 
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*(•>»'*) - X I A)m 6in jf 


k=l n=l 


e ikx e O)t 


(M) 


and producing calculations analogous to those which were mentioned above 
for the general case n ^ 0, gives the formula for the critical velocity 
of flow, as follows: 



Formula (4.5) coincides with the critical velocity of flow for an 
infinite plate, simply supported along the edge parallel to the direc- 
tion of flow, which was obtained first in reference 4 . 

^Omitting all intermediate calculations, we quote the formula of 
the critical velocity of flow for the unbounded closed cylindrical shell, 
found from consideration of the differential equation of small vibrations 
of cylindrical shells of medium length (2.4): 



It is possible to make use of formula (4.6) for all values of 
n ^ 2 . From this formula it follows, that for n = <», the critical 
velocity of an unbounded closed cylindrical shell of medium length coin- 
cides with the velocity of the unbounded closed cylindrical shell found 
starting from the theory of shallow shells . 


^Translator ’ s note: This paragraph and. the one following seem mis- 

placed. They were probably intended to conclude the preceding section. 
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5. FLUTTER OF CLOSED CYLINDRICAL SHELLS OF BOUNDED LENGTH 
FOR DIFFERENT BOUNDARY CONDITIONS ON THE ENDS 


The equation of small vibrations of shells of medium, length (2.4) 
is used to examine a series of boundary problems . In this examination 
is introduced a new variable 6, connected with a by the formula 

a- = I (5-1) 

Then in a form convenient for solution, the equation of small vibration 
(2.4) is written 

^Vl + c 2 iVd 2 1 l\ 2 ^*1 1 p ii_ _^fl_ + B^_ d^ $ l bvz 3 d^q _ 0 
^ * R^Vdp 2 / dp 1 * - ER 2 dt 2 ^ EhR 2 dt dp^ EhR 2 86 dp 1 " 

(5.2) 

To equation (5*2), in each particular case, must be added boundary 
■conditions on the ends 1=0 and 6 = 1. 


By defining, according to formulas (1.6) and ( 1 . 7 ), the displacements 
and internal forces of the shell through the boundary conditions for 

the boundary problems may be represented in the following forms: 


(a) For a shell simply supported on the ends | = 0,1 = 1: 


w = 


th., 0 

dp 1 " 


Mq = £ cr 


8^ $2. <^$ 2 . 


R L8p 6 8p 4 J 

(b) For a shell clamped on the ends 1=0 and 1=1: 


=0 at |=0 and £ = 1 

(5-3) 




w = 


1 = 0 dw = R ^<Eq = 0 


ar 


d6 1 


86 dp" 


at 6=0 and 6=1 

(5-4) 


(c) For a shell simply supported on the end 6=1 and rigidly 
clamped on the edge 6=0: 
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w = 


= 0 

dw 

_ R 

d 5 tl 

at 

S = 0 

"s 






as 

Z 

as ap 4 - 




(5.5) 











W a 

^41 

= 0 


_ D 

a 

d 6 $>i a^«i 

= 0 

at 

s = l 




dp 1 *" 


R 


ap 6 ap 4 _ 



> 




(d) 

For a 

shell clamped on the end i = 

0 and 

force on the 

end 

1 = 

1: 











W = 


= 0 

dw 

_ R 

m o 

at 

S = 0 

\ 






as 

' Z 

as ap 4 




(5-6) 











N, = 


= 0 


S 

_ EhR 2 o 4- *! 

- = 0 

at 

S = 1 



a. 

Z 2 

as 2 ap 2 




z 5 as 5 ap 


J 




(From, the second group of relations it is seen that the "boundary con- 
ditions on the free end are partially satisfied . ) 


(e) For a shell simply supported on the end 1=0 and free on 
the edge 1=1: 


w = 


d 4 «h D 

a 6 <ti 8^4>i 
■ + , 

1 E 

Lap 6 ap 4 _ 

_ HUE a 4 *! . 0 s 

EhR 2 ^ $ 1 


l 2 d£ 2 d3 2 


= 0 


= 0 


at i = 0 


at | 


( 5 . 7 ) 


For the class of solutions 

00 00 

£ CknX k (|) e^cos n3 t (5*8) 

n=2 k=l 




equation (5*2) is written after a series of simple transformations, 
in the form 


d 4 X k k 
— - An 4 


L2( n 2 

d| L 1 v 


l) n 4 - An 4 ! X k = 0 


4 


( 5 . 9 ) 



NACA TM 1438 


15 


where 



(5-10) 


The equation of the parabola of stability will have the form: 


X 1 = p 


h 2 ER 2 ~y 2 

~2lT 

B f-l* 




( 5 . 11 ) 


For fixed c-j_, n, A, and \ the solution of equation ( 5-9 ) , in 

the case when the roots of the characteristic equation are distinct, has 
the form 


X k (|) . C ie - k ^ + C 2 e _k21 + C 5 e" k51 + C^e"^ (5-12) 


The rest of the problem reduces to the determination of nontrivial 
solutions C± } for which it. is sufficient to subject the solution 
(eq.. ( 5 . 12 )) to the boundary conditions and to require the vanishing 
of the appropriate determinant A(kj_ ) . Avoiding the problem of the 

form of the determinant A(kj_) for the different possible combinations 
of multiple roots, the function 


r(ki) 


'6(ki) 


( 5 . 13 ) 


is introduced where 

8 ( k l) s ( k l “ k 2)( k l - k 3)( k l - *4)0*2 - k 3)( k 2 - ( k 3 - * 4 ) 

From the expression 6 (kj_) it follows that all zeros of the func- 
tions A^k^) will exist by virtue of zeros of &(kj_), and F(kj_) will 

be an analytic function in the whole domain of the variation of the 
variables . 
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The solution of equation (5*9) in its general case is accompanied 
with considerable mathematical difficulties, and for this reason we apply 
here the method of investigation of the eigenvalues A which was proposed 
in references 4 and 5. The essence of the method is contained in that, 
instead of the solution of equation (5.9) , the parameters of the prob- 
lem A, A and the two sought for roots k^,k^ are expressed through 
two other roots k^_, kg of the equations — 


A = ■ + *i2 2k l + k i 2k 2 + k J\ 

A = c-j_ 2 (n 2 - l) 




,2 kikgkjlq,. 


(5.14) 


k 3,b = s ~ * (“A - $(h * ^ 2 )^ 


1/2 




and instead of finding the eigenvalues of equation (5*9) there is investi- 
gated a system of two equations comprising the characteristic system 


A + - 7*)= 0 

r 

F (q, 7 ) = t^hlL =, 0 
s(n , 7 ) 

where r\ and 7 are quantities connected with the roots 


(5.15) 


k]_ = ti + ly kg = T) - iy (5.16) 

of the equation and 

8 (t), 7 ) = l 6 i 7 [t 2 - 2tj 2 J ^ [(t 2 - 3 t] 2 ) + 4 ti 2 7 2 J (5.17) 

The left part of each of equations (5-15) represents, an analytic 
function of the variables tj arid 7, and the problem consists in finding 
a solution such as 

T)1 = Tlt(n,A) 7± = 7i( n ,A) (5.18) 


of the system, which permits, through use of the formulas 
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rr" 

k 3> - -Ilf} 2 - 2 1 ! 2 ] ^ 

* = ci 2 (i, 2 - l) 2 + 3 !_l£(,2 . 3, 2 ) 

nr 


(5.19) 


for each boundary problem, to compute the corresponding eigenvalues A 
and to establish that value A at which an eigenvalue becomes complex. 

The solution of the characteristic system is always easier to obtain 
graphically if graphs of the curves defining the equations (5-15) are 
penciled on one sheet in a rectangular system of coordinates q, 7. The 
general form of the curves of the characteristic system is shown in fig- 
ure 2; the graphs of the curves corresponding to the first equation of 
the system (hyperbolas) axe drawn at different values A =» Constant. 

The rest of the problem reduces to the establishment of those values 
A*i where points of the first and second real branches (5.18) coincide, 

'and it is impossible to make any deduction concerning the eigenvalues 
of the examined boundary problems . 

Equating A = A^ according to equations (5. 10 ), the velocity of 

flow at which there still exists stability of the unperturbed motion, 
but above which the motion may possibly become unstable is found. Con- 
sequently, for each particular boundary problem it is necessary, first 
of all, to compose the expression for the second equation of the charac- 
teristic system a(t\, 7) = 0 . 

The composition of the determinant A(t\, 7) is shown for the example 
of a simply supported shell. For the determination of nonzero 
C-^(i = 1, 2, 3, 4 ), expression (5*12) for Xj^Ct) is subjected to the 

boundary conditions (5*3) and the determinant of the system thus obtained 
is equated to zero: 


1 





= 0 






k 3 2 e' k 5 
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The determinant having been disclosed, producing in it a replacement 
of in terms of ti and 7 according to the formulas (5.16) results 

in 


(a) 


^,7) = i 


-2r\ 2 yfj 2 - 2r\^j ^ cosh 2ti + 2t\ 2 7 [j 2 - 2t^} ^ cos T^oshjj 2 - 2 tj^} ^ 


(3^- 7^- 2t| 2 7 2 ) sin 7 sinh fy 2 - 2 t\ 2 \ ^ ll6i = 0 


( 5 - 20 ) 


) 


The expression A(t|, 7) for different boundary problems is obtained 
by analogous means: 

(b) For a shell clamped on the ends £ = 0 and £ = 1, 

A(tj, 7) = 8i< 7 [t 2 - 2tj 2 J ^ cos 7 cosh^T^ - 2 ti 2 ) -cosh 2 t)J + 


3ri 2 sin 7 sinhjT 2 - 2 ti 2 J j* = 0 


(5.21) 


(c) For a shell clamped on the end £ » 0 and simply supported on 
the end £ = 1, 


4(i|,7) = 8i 


2 ti 7 It 2 - 2 t} 2 J sinh 2r\ + 


(t 2 ~ Jf) 2 ) (t 2 - 2r) 2 J 1 ^ 2 sin 7 cosh {t 2 - 2i]^\ ^ - 


y(j2 + ri 2 ) cos 7 sinh [t 2 - 2ti 2 | 


1/2 


= 0 


( 5 - 22 ) 


(d) For a shell damped on the end £ = 0 and free on the end £ - 1, 
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A(t], 7) = ij^Ol 2 + 7 s ) 2 [ t 2 ~ 2ti 2 ] ^cosh 2 t \ + 

1^7(26^ + 27^ - 4t| 2 7 2 )[7 2 - 2ti 2 3 1 / 2 cob 7. cosh, [t 2 - 2tj 2 ] + 

3-/2 

8t| 2 (2ti 2 t 2 - 7^ +, 3^) sin 7 sInhjT 2 - 2 t} 2 ] 

16^7(7^ - t] 4) cos 7 sinhjVs - 2^ ' - 

L/2 p — I !L/2 

l6r] (l+^T 2 - 3^ - 7*0 [j 2 - 2 t) 2] sin 7 cosh [t 2 - 2t] 2 J 

32 t) 2 7 2 (t 2 - t\ 2 ) [t 2 - 2t\^] e" 2T ^| = ® (5-23) 

(e) For a shell simply supported on the end £ = 0 and free on 
the end £ = 1, 

A(H,7) = i|- 2 ri 7 ( 7 ^ + t^) [ r 2 _ 2 n + 

T}7 (t 2 - 2^2 ) ^ j/y 2 - T] 2 ^ + (72 _ 3^2) J e“ 2T ^ + 

81^7 [t 2 - 2^ 1//2 (7 2 - n 2 ) cos 7 cosh [y 2 - 2 t^] ^ + 

^n(3n 2 7^ - 76 + 3n 6 - ^T 2 ) sin 7 sinhlT 2 - 2r\^ ^ + 

7 { 5 ^ 2 7 ^ - 7 ^ - 19V* - T 2 + 23 ^ cos 7 sinhjjy 2 - 2r\2j ^ + 

[72 - 2rj2 ] 1//2 (76 + Ht i^t 2 - t^ 2 7^ - 3tl 6 )sin 7 cosh [t 2 - 2rfi\ = 0 

(5-24) 
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It is noted that, for rj = 0, the equation A(r), 7) = 0 degenerates 
into the characteristic equations of the team functions fundamental for 
the corresponding boundary problems. 

For all the considered boundary problems, the graphs of equations 
(5-15) were constructed. With the help of these, the critical velocities 
of flow (fig. 2), that is, the velocities above which ..the unperturbed 
motion of the shell becomes unstable, were determined. Values of the 
critical velocities for three boundary problems (shells, simply supported 
on two sides, clamped on two sides, and clamped on one_ and simply sup- 
ported on the other side) are quoted in table 1, for different ratios 
R/Z and h/R. 

Computation of the critical velocities for the cases shown were 
conducted for n = 4 at the following values of the constants: 


E = 2 X 10 12 P a 7.8 0 = 0.3 

cm? cm? 

p 0 = 1.014 x 10 ^ v 0 a 3.4 x 10 ^ k = 1.4 

cm? 

Here will be considered first the case where, for the class of thin 
shells of short and medium length, four transverse half-waves correspond 
to the fundamental mode of free vibrations of the one longitudinal half- 
wave . (The circle passes to an ellipse . ) 

It is possible to establish this situation if the frequency of free 
vibrations of closed cylindrical shells for the case of hinged support 
on the edges is investigated. 

The minimum frequency of free vibrations of a simply supported 
closed cylindrical shell, calculated according to the theory of shallow 
shells, occurs at conditions coinciding with conditions (3.21), at which 
occurs the minimum velocity of propagation of a traveling wave along the 
generators of the shell in the absence of flow. 

Because one of the equations of the characteristic system F(t],7) = 0 
does not depend on n, then from the first of equations (5*18) it follows 
on first glance that, for other conditions equal, with increase of the 
number of terms of the expansion n, self-induced vibration can possibly 
occur at lower velocities of flow. It is seen that there exist waves, 
the velocity of propagation of which (for all n, larger or smaller than 
n = 4 ) is above the velocity of propagation of the traveling wave occur- 
ring for n = 4 . This situation has been successfully corroborated only 
for the case where the velocity of flow equals zero. 
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6. APPLICATION OF THE METHOD OF EUBNOV-GALERKEN TO THE 
INVESTIGATION OF SELF- INDUCED VIBRATIONS 
OF CYLINDRICAL SHELLS 


In the preceeding paragraph a method was presented for the investi- 
gation of the eigenvalues of boundary problems resulting from, the equa- 
tion of small vibrations of cylindrical shells of medium length, which 
permits the determination of exact values of the critical velocities of 
flow. 


An analogous method could not be successfully applied to the invest! 
gation of the eigenvalues of the separate boundary problems using the 
general equation of small vibrations of cylindrical shells (2.3). There- 
fore, for consideration of problems on self-induced-vibrations of closed 
cylindrical shells and cylindrical panels according to the theory of 
curved shells, a variational method was applied. First, there were deter 
mined by the variational method, the critical velocities of flow for the 
class of closed cylindrical shells of medium length with the different 
boundary conditions on the ends. 

The values of the critical velocities resulting from the equation 
of shells of medium length, having been found in the second and third 
approximations, are given in table 1. From the table it is seen that 
the second approximation according to Galerkin gives a somewhat lower 
value of the critical velocity and the third approximation a somewhat 
higher value of the velocity as compared with the exact value of V^j 

that is, the second and third approximations bracket the exact value 
of the critical velocity. 

The good convergence of the variational method permitted its appli- 
cation to the investigation of problems of the self-induced vibration 
of closed cylindrical shells and cylindrical panels using the more 
general equation of shallow shells (eq. (2.3))* 

In table 2 are quoted values of the critical velocities for a 
closed cylindrical shell simply supported on the ends, having been 
found starting from the theory of curved shells . The information given 
in table 2 with the corresponding magnitudes of the critical velocities 
quoted in table 1 makes it possible to remark that the difference in 
magnitude of the critical velocities does not exceed 10 to 15 percent 
and, consequently, for practical purposes it is entirely Justified for 
the investigation of problems of the self-induced vibration of closed 
cylindrical shells to make use of the simple equation of shells of 
medium length. 
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TABLE I— VALUE OF THE CRITICAL VELOCITIES (m/sec) OF 
FLOW FOR CLOSED CYLINDRICAL SHELLS 

[shells simply supported on the ends ( I) ; clamped on one 
end and simply supported on the other (ll)j 
and clamped on "both ends (ill)] 




Exact 

By the variational method 

h 

R 

R 

l 

T 

II 

III 

2nd approximation 

3rd approximation 



X 

I 

II 

in 

I 

II 

III 


1 

6 

lit-, 067. 

19,182 

26,642 

12,278 

16,883 

22,132 

14,770 

20,448 

29,040 

1 

1 

8 

5,934 

9,092 

11,239 

5,203 

7,141 

9,352 

6,332 

8,869 

11,880 

200 

JL 

10 

3,038 

4,143 

5,754 

2,708 

3,686 

4,812 

3,369 

4,238 

6,134 


1 

12 

1,758 

2,397 

3,330 

1,648 

2,186 

2,822 

1,999 

2,656 

3,526 


1 

6 

9,378 

12,788 

17,761 

8,197 

11,266 

14,765 

9,847 

13,939 

19,360 

1 

1 

8 

3,956 

5,396 

7,493 

3,484 

4,775 

6,249 

4,308 

5,913 

7)920 

300 

1 

10 

2,025 

2,769 

3,836 

1,687 

2,475 

3,225 

2,246 

3,314 

4,089 


1 

12 

1>472 

1,598 

2,220 

1,120 

1,444 

1,902 

1,390 

2,014 

2,539 


1 

6 

— 

7,033 

9,591 

13,321 

6,160 

8,46l 

11,085 

7,385 

10,723 

14,879 

\ 

i 

2,968 

It-, Oit-8 

5,623 

2,629 

3,596 

4,701 

3,233 

4,437 

5,943 

boo 

i. 

10 

1,520 

2,072 

2,879 

. 1,378 

1,874 

2,436 

1,685 

2,533 

3,069 


1 

12 

879 

1,199 

1, 666 

861 

1,130 

1,447’ 

1,043 

1,626 

1,976 


1 

6 

5,626 

7,673 

10,656 

4,940 

6,781 

8, 880 

5,908 

1 

8,578 

| 

H,903 

1 

500 

1 

5 

2,375 

3,238 

4,498 

* 

2,119 

2,892 

3,776 

2,586 

3,549 

4,754 

1 

10 

1,216| 

1 

1,658 

2,303 

1,129 

1,512 

1,968 

1,348' 

2,026 

2,595 


1 

z 

3,751 

5,115' 

7,104 

3,322 

4,548 

5,946 

3,938 

5,718 

7,367 

1 

750 

1 1 
8 

1,583 

f 

2,159 

2,999 

1,450 

1,964 

2,552 

1,890 

2,366 

3,169 

1 

10 

810 

1,105 

1,535 

796 

1,056 

1,354 

1,025 

1,350 

1,636 




h 

IT 


E 

l 


1 

200 


1 

300 


1 

400 


1 

500 


1 

750 


By 2nd approximation 
m/sec 


10, 858 
4,959 


_ 1 _ 

10 

JL_ 

12 


1 

6 

1 

8 

_ 1 _ 

10 

_1 

12 


1 

6 

l 

8 

_1_ 

10 

1 

6 

1 

8 

j__ 

10 

12 


1 

6 

1 

8 


2,684 

1,689 

7,185 

5,261 

1,770 

1,117 

5,584 

2,449 

1,557 

4,315 

1,971 

1,086 

700 

2,900 

1,548 


By 3rd approximation 
m/sec 

14,770 

6,498 

5,463 

2,235 

9,847 

4,352 

2,208 

1,452 

7,585 

5,167 

1,752 

5,908 

2,456 

1,268 

956 

5,958 

1,775 


_l_ 

10 


765 


1,050 
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